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ABSTRACT 

An attempt has been made to model the so called “Leading Edge Vortex” 'which 
forms over the leading edge of delta wings at high angles of attack. 

A simplified model has been considered, namely that of a two-dimensional, 
inviscid, incompressible steady flow around a flat plate at an angle of attack with a 
stationary vortex detached on top, as well as a sink to simulate the strong spanwise 
flow. 

The results appear to agree qualitatively with experiments. A comparison has 
also been made between the lift and the drag of this model and the corresponding 
results for two classical solutions: 

(i) that of totally attached flow over the plate with the Kutta condition satisfied at 
the trailing edge only, 

(ii) the Helmholtz solution of totally separated flow over the plate. 
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NOMENCLATURE 


English 

a 

b 

e 

C D 

Cl 

Cm 

Cp 

d 

D 
F 
k 
k o 
k i 
K 
K o 
K x 
L 

rn 
m i 
Mi 

n 

r,0 
R i 

S 

u,v 

Voo 

w 

x,y 

z 

Z 

z' 


letter symbols: 

radius of the cylinder, 
length of the span of the plate, 
chord length of the plate, 
drag coefficient, 
lift coefficient. 

pitching moment coefficient, 
pressure coefficient. 

distance between the streamlines leading to the two stagnation points on 
the pl^te. 
drag force. 

total force on the plate. 

total vortex strength. 

bound vortex strength. 

leading edge vortex strength. 

non-dimensional total vortex strength. 

non-dimensional bound vortex strength. 

non-dimensional leading edge vortex strength. 

lift force. 

mass flow rate. 

sink strength. 

non-dimensional sink strength. 

non-dimensional parameter showing position along the plate, 
porlar coordinates in the z-plane. 

non-dimensional distance from the center of the cylinder to equilibrium 

point. 

wing area. 

velocity components, 
free stream velocity, 
complex potential. 

cartesian coordinates in the 2- plane, 
complex variable in the original circle plane, 
complex conjugate of z. 
complex variable in the rotated circle plane. 


Greek letter symbols: 
a angle of attack. 

T total circulation. 

To circulation due to the bound vortex. 

{ complex variable in the Gnal plate plane. 

{' complex variable in the Joukowski plane, 

rj coordinates in the {-plane. 

r\' coordinates in the {'-plane. 

Ao doublet strength. 

p air density. 


Subscripts: 

c in the plane of the cylinder. 

p in the plane of the plate, 

r radial component. 

sc on the surface of the cylinder. 

sp on the surface of the plate. 

TE at the trailing edge. 

z in the 2 -plane. 

6 tangential component. 

{ in the {-plane. 

0 refers to the origin (center of the cylinder). 

1 refers to the equilibrium point. 
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1* INTRODUCTION 


So far there are two well-known models for the flow over a flat plate at an angle 
of attack. That of totally attached flow with the Kutta condition satisfied at the trailing 
edge only, and that of totally separated flow (Helmholtz solution). 

The present model, considering partially separated flow lies somewhere between the 
two and despite the fact that it too is two-dimensional, gives some very useful represen- 
tation of the three-dimensional flow over delta wings. On such wings the leading edge is 
usually sharp, causing thus flow separations even at moderate angles of attack. These flow 
separations take the form of two free vortex layers joined to the leading edge of the wing 
and rolling up to form spiral shaped vortices above the upper surface of the wing (Figure 
la). 

These vortices induce additional velocities at the upper stirface of the wing. The 
corresponding pressure distribution shows distinctly marked minima beneath the vortex 
axes (Figure lb). Accordingly, an additional lift force occurs which depends nonlinearly 
on the angle of attack (Figure lc). 

In the theoretical study made here, a simplified model has been considered, namely 
that of a two-dimensional, inviscid, incompressible steady flow over a flat plate at an angle 
of attack, with a stationary vortex detached on top ns well as a sink to simulate the strong 
spanwisc flow caused by the pressure gradient due to sweep in the three-dimensional case. 
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2. ANALYSIS OF THE FLOW FIELD 


In the original circle plane (also called 2 -plane) the How field consists of the following 
components: 

(i) Uniform wind coming from the negative z-axis. 

(ii) A doublet A<> at the origin to simulate a circular cylinder |z| ss a, 

(iii) A bound vortex ko at the origin to account for the circulation To around the plate. 
Note that although the flow is aligned with the z-axis the plate is at an angle of 
attack, requiring thus circulation To for the Kutta condition to be satisfied (sec 
Figure 2a). 

(iv) A steady vortex k\ of finite radius placed at zi(ri,0i) to simulate the lending edge 
vortex. 

(v) A steady vortex — fci at thf diverse square point of Zi induced by the circle theorem 
(see Appendix Al.I). 

(vi) A steady vortex ki at the origin also induced by the circle theorem (sec Appendix 
Al.I). 

(vii) A sink mi (mi < 0) placed at zi (n,0i) to simulate the spanwisc flow along the 
center of the vortex. 

(viii) A sink mi at the inverse square point of zi, induced by the circle theorem (see 
Appendix A1.2). 

(ix) A source —mi at the origin also induced by the circle theorem (sec Appendix Ai.2). 
2.1 The Complex Potential 

For the components described above (regrouping them together), the complex po- 
tential is given below: 


V) 


= Voo +(“ m i+*fc) ln#+(mi+sfci) ln(z- 2 i)-f (mi-t'fci) In ^ 2 - (1) 


where 

k — fco •+• ki (2) 


2 


2.2 Th« velocity Field 

Differentiation of * he convH'X potential gives the velocity field: 

~ a (u r ~iuo)e~ i0 (3) 

which applied to equation (1) gives: 


2.2.1 Radial Components 


«r 


f oo 


+ 



miri[rn 


r 


a m i , mi[r - ri cos (0 - 0i)] + kin sin(fl - 0i) 
C0S r r 2 + r 2 - 2rri cos(0 — 0|) 

— a 2 cos(0 - 0|)| — fc ^a 2 sin(fl — 0x) 

2 rf -f a 4 — 2 tti<i 2 cos(0 —Oi) 


( 4 ) 


2.2.2 Tangential Components 


Me 


r OO 


+ “2 ) 


A m ifi sin(0 - Oi ) - fci[r ~ n eos(fl - flj)j 
r r 2 + rf ~ 2rrji coa(0 — 0i) 


+ 


mino 2 sin(fl — 0i) 4- fejri[rri — a 2 cos(fl — fl|)] 
r 2 rj -fa 4 — 2rrja 2 cos(0 — 0|) 


( 5 ) 


2.3 Conformal Transformation Used to Analyze the Flow Field 

A solution for the flow field is provided by a mapping sequence that transforms the 
potential flow about a circle into a flow about a flat plate at an angle of attack with a 
detached vortex/sink combination. The steps are the following (sec also Figure 2): 

a. original circle plane (z-plane), 

b. rotated circle plane (z'-planc), z' — ze ,a , 

c. Joukowski plane (f'-planc), = z' + fr» 

d. final plate plane (^-plane), = f'e”'®. 

The first step is a simple rotation in order to make the flat plate depicted inside the circle 
aligned with the z-axis. The second step is the Joukowski transformation which transforms 
the circle into a flat plate. The third step is another rotation which gives the plate an 
angle of attack with respect to the horizontal free stream velocity. 


( 6 ) 





Combining tlio three steps wo have that 


x 


or 


and if we set 


r=-± 

2 


^(|) a -(«■ 


•»oi 


f = £ 4- «i? 


we get 


£ = r cos 0 4 cos(0 4* 2o) 

r 

a 3 

q = r sin 0 — — - sin(0 4- 2a) 


On the surface of the cylinder (corresponding to the surface of the plate) r — a so 

£«p — a [cos 0 4* cos (0 + 2a)] 

W»p = « [sin 0 — sin (0 + 2a)] 


now 


and 


dw _ dw dz 
d$ dz d$ 


( 6 ) 

( 7 ) 

( 8 ) 

( 9 ) 

( 10 ) 

( 11 ) 

( 12 ) 

(13) 


dz 


d$ d$/dz 


z 3 e 3ia 


z 2 e 2ia — o 2 
r 2 t i(2a+20) 


- r 2 e i{2a+20) _ 0 2 

Considering the surface again (r = a), expanding the right-hand side into sines and cosines 
and using equation (3) we get finally 

r -I _ sin (40 + 4a) — 2 sin (20 4- 2a) — 2 cos (20 4- 2a) , . . . 

'■"'••'U 2 4-2 cos (40 4- 4a) - 4 cos (20 4- 2a) 1 ' 

1 




2 cos (20 4- 2a) 


KJ, 
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(14) 

(15) 


3. CONDITIONS 


Boundary conditions in the final plate plane require that: 

(i) the flow depart smoothly from the trailing edge; 

(ii) the flow depart smoothly from the leading edge, i.e., the Kutta condition must 
be satisfied at both edges of the plate. 

Also, the vortex/sink combination must be located at an equilibrium point in the 
flow field, that is, a point where the velocity induced by all other components is zero. 

3.1 Kutta Condition at the Trailing Edge 

Prom equations (14) and (15) with 0 = -dt we find the velocity components at the 
trailing edge of the plate: 

kj” = -J ( 16 ) 

K.]”- < 17 > 

For the Kutta condition to be satisfied there (i.e., u r finite and uo = 0) we must have 

[« 0 jr -a =o (is) 

or, from equation (5) with r = a and 0 — —a 
(‘iVooOsina — k) [r 2 + a 2 — 2na cos(0i -f a)] — 2mirjasin(0i + a) + &i(r 2 — a 2 ) = 0 (19) 

3.2 Kutta Condition at the Leading Edge 
Similarly, for the leading edge we must have 

K.1T- -• P») 

and again from equation (5) with r — a and 0 — re — a 
(2 Food sin a + k) [r 2 -f a 2 + 2riacos(0i -|-a)] - 2miria sin(0i -t-a) - fci(r 2 - a 2 ) = 0 (21) 
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3.3 Vortex/Sink Velocity Condition 

Equation (3) cannot be used to find the conditions which make the velocity vanish 
at the location of the sink and vortex because the velocity given by equation (3) is singular 
there. Threfore, the usual limiting process has to be used to find the velocity components 
at the center of the vortex/sink combination. The strengths of the singularities can then 
be adjusted so that they arc stationary in the presence of the plate. Following the analysis 
of reference 3 the complex velocity function in the final plate plane at the equilibrium 
point is found to be 


[u p 



1 

2 


(mi + iki) 


[dz/d<;] {l 


= 0 


The derivative ^ has been found in Section 2.3: 

dz _ (»e ia ) 2 

d$ (ze'“) 2 — d~ 


( 22 ) 


while 

d?z _ d / dz\ _ d T z 2 e 2,a 1 
d? 2 df \df/ d$ [z 2 e 2,a — a 2 J 

and substituting z from equation (7), taking the derivative and substituting back £ from 
equation (6) we get 

d 2 z _ a 2 e 2,, *[4 z A + 3(z 2 -j- a 2 e -2, “) 2 j 
df 2 2 z (a 2 e -2, “ — z 2 )(z 2 e 2, “ — a 2 ) 2 

Also, from Figure 3 we have the following transformations: 


ti r = u c cos 0 + v c sin 9 
uq = v c cos 0 — u c sin 0 


or 

. u c = u r cos 0 — uo sin 0 

v c = u r sin 9 -(- uo cos 9 

Substituting all these into equation (22) we get 


[(u rl cos 9i - uo i sin<?i) - * (u r i sin0i + uoi cos0i)] z 2 e 
/_ | -*jl X H + 3 (*i + a 2 e -2, “) 2 ] _ 




( 23 ) 


where 


zt — rie ,tf * = ri (cos 0i + i sin 0i) (24) 

Splitting into real and imaginary parts and substituting for it r i and uqi their equivalent 
expressions from equations (4) and (5) with r = n and 0 = 0j we get 

21^ r i(rJ — a 2 ) 2 cos0j.cos(60i 4* 2a) — 4V 00 rJa 2 (rJ — a 2 ) 2 cos0x cos40i 

+ 2Vooria 4 (r 2 - a 2 ) 2 cos0i cos(20i - 2a) - 2V’ 00 r®(r 4 — a 4 )sin0i sin(60i •+• 2a) 

+ 4V’oo^i<* 2 (»’i — ° 4 ) s ‘ n s ‘ n 401 — ZVoofia 4 ^ 4 — a 4 ) sin 0i sin(20i — 2a) 

- 2 fc 0 r®(r 2 - a 2 ) sin(60i + 2a) + 4&or 4 a 2 (r 2 - a 2 ) sin 40i 

- 2fcor 2 a 4 (r 2 - a 2 ) sin(20i - 2a) + 2 &ir®a 2 sin(60i + 2a) 

+ &ir 4 a 2 (7r 2 — 11a 2 ) sin 40 1 + 7Aqr 4 a 2 (r 2 — a 2 ) cos 40i 

+ 2fcir 2 a 4 (3r 2 - 2a 2 ) sin(20i - 2a) + 6A:ir 2 a 4 (r 2 - a 2 ) cos(20i - 2a) 

+ 3(ki — mi)(rj — a 2 )(cos 4a — sin 4a) a 0 -f- 2mir®a 2 cos(60i + 2a) 

- mir 4 a 2 (7r 2 — 3a 2 ) cos 40i - 7mir 4 a 2 (r 2 — a 2 ) sin 40i 

- 2mir 2 a 4 (3r 2 - 4a 2 ) cos(20i - 2a) - 6mir 2 a 4 (r 2 - a 2 ) sin(20i - 2a) = 0 (25) 

and 

2V’ 00 »’i(r 2 — a 2 ) 2 cos 0i sin(60j + 2a) — 4V 00 r|a 2 (rJ — a 2 ) 2 cos 0i sin 40 1 

+ 2Vooria 4 (rJ - a 2 ) 2 cos 0 X sin(20i - 2a) + 2V' 00 rJ(rJ - a 4 ) sin O t cos((i0 t f 2a) 

- 4V’ <x >r 2 a 2 (r 4 — a 4 ) sin0i cos40x + 2V’ 00 ria 4 (r 4 — a 4 )sin0i cos(20i ~ 2a) 

+ 2&or®(r 2 — a 2 ) cos(60i + 2a) — 4fc(>r 4 a 2 (r 2 - a 2 ) cos40i 

+ 2A: 0 r 2 a 4 (r 2 - a 2 ) cos(20i - 2a) - 2fci r®a 2 cos(60i + 2a) 

— A?ir 4 a 2 (7r 2 — lla 2 )cos40i — 7&ir 4 a 2 (r 2 — a 2 ) sin 40* 

- 2A:ir 2 a 4 (3r 2 - 2a 2 ) cos(20x - 2a) - dkirja 4 (rj - a 2 ) sin(20i - 2a) 

- 3(&i 4- mx)(rj — a 2 )a°(cos 4a — sin 4a) -f 2mir®a 2 sin(60i -f- 2a) 

- mir 4 a 2 (7r 2 - 3a 2 ) sin 40i - 7mir 4 a 2 (r 2 - a 2 ) cos 40i 

- 2mir 2 a 4 (3r 2 - 4a 2 ) sin(20i - 2a) - 6mir?a 4 (rjf - a 2 ) cos(20i - 2a) = 0 (26) 
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3.4 Non-Dimensional Parameters 

At this point some non-dimensional parameters have to be introduced if the system 
of equations (19), (21), (25), and (26) is to be solved explicitly. These are the following: 


If - k ° 

Ko ~^z 

(27) 

K kl 

Kl 

(28) 

*■=-£- 

dVgg 

(29) 

iu - mi 

(30) 

R t = — . 
a 

(31) 


Substitution into the original equations yields the following system of four equations 
in five unknowns Ko, ffi, Mi, ^it and 0\\ 

(2 sin a - K)[R\ + 1 - 2 Ri cos(0i + a)] - 2MiR\ sin(0i + a) + K^R\ - 1) = 0 (32) 
(2 sin a + K)[R\ + 1 + 2f?i cos(0 A + a)] — 2 Mi Ri sin(0i + o) — Kx(R? — 1) = 0 (33) 

2/2i(f?i — l) 2 cos 0i cos(60i -|- 2a) — 4 Ri(R? ~ l) 2 cos0i cos40i 

+ 2Ri{R\ - l) 2 cos 0i cos(20i - 2a) - 2R\{R\ - 1) sin 0i sin(60i + 2a) 

+ iR^R 4 — 1) sin 0i sin 40i - 2Ri(R 4 — 1) sin0i sin(20i — 2a) 

- 2 Jfoitf (A? - 1) sin (60 1 + 2a) + 4 K {) R\{R\ - 1) sin 40i 

- 2tf 0 fii (^i - 1) sin (20i - 2a) + 2K l R% sin(60i + 2a) 

+ KlR\ {1R\ - 11) sin 40i + (Rf - 1) cos 40i 

+ 2KiRl{ZRl - 2) sin (20i - 2a) + dKiR]{R\ - 1) cos(20i - 2a) 

+ 3(Ki - Mi)(f? 2 — 1) (cos 4a — sin 4a) 4 2Afi/?i cos(60i + 2a) 

- MiR^lRl - 3) cos 40i - 1MiR\{R\ - 1) sin 40i 

- 2M l R\{ZR 2 i - 4) cos(20i - 2a) - 6Mi/i?(fi? - 1) sin(20i - 2a) = 0 (34) 

2R\{R\ - l) 2 cos 0i sin(60i + 2a) - 4ft?(7 1\ - I) 2 cos0isin40i 

+ 2Ri{R\ — l) 2 cos 0i sin(20i — 2a) + 2R\(R 4 — 1) sin0i cos(60i 4- 2a) 
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- 47?J(i2| — 1) sin 0i cos 40i + 2 Ri(R* — 1) sin 0i cos(20i — 2a) 

+ 2KoR°(Rl ~ 1) cos(60i + 2a) - AK^R^Rl - 1) cos40i 

+ 2K{)Rl{R\ - 1) cos(20i - 2a) — 2Kii2® c os(60i -f 2a) 

- KiR\{lR 2 x - 11) cos 40! - IKiR^R] - 1) sin 40! 

- 2K i R\{ZR\ - 2) cos(20i - 2a) - 6KiRj{Rj - 1) sin(20i - 2a) 

- 3(Ki -H Mi)(R* — l)(cos 4a -- sin 4a) + 2 MiR* sin(60i + 2a) 

- M t Rj(7R 2 L - 3) sin 40 x - 7 M x R\{Rl - 1) cos40j 

- 2MiR 2 (ZR J — 4) sin (20i — 2a) — 6MiR 2 (R 2 — 1) cos (20 1 — 2a) = 0 (35) 


3.5 Additional Condition 

We see that the equations derived so far are not enough to give us a unique solution. 
Therefore we must seek additional information in the nature of the flow. As mentioned 
in the introduction the flow is assumed inviscid which implies that the total force on the 
plate must be perpendicular to it (sec Figure 4). Therefore 

tan a = — (36) 

Lt 

Now we have to relate the drag <uid the lift with the unknowns K and Afi. 

From Figure 5 we sec that the fluid between streamlines a and b dissappears into the 
sink producing drag which can be calculated ns follows: The rate of mass of fluid between 
the streamlines a and b (which is distance d apart) must be equal to the rate with which 
mass is dissappearing in the sink, i.e., 


m = -pV^db — mipb2n. 


(37) 


Now the drag is given by 'the rate of loss of momentum of that fluid therefore 


D = -mVoo 


combining equations (37) and (38) we get 


mi 


D 


2npbV ( 


oo 


(38) 


(39) 
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and from equation (30) we get 


(40) 


Mi = 


D 

2npbaV£, 


but 


S = be 

and from the Joukowski transformation we know that 


c = 4o 


also, 


Cd — 7 


D 


hpVlS' 


The last four equations combine to give 


C D = -ttMi 


The lift can be found from the Ku tt a- J oukowski law 


and since 


and 


we get 


L = pV ooT • b 

r = 2wk 
L 


C i = T 


kpVis 


c l = *k 

alter employing equation (29). 

Now going back to equation (36) we sec that 

Mi 


tana = 


K 


(41) 

(«) 

(43) 

(44) 

(45) 

(40) 

(47) 

(48) 


(49) 


thus we have a system of five equations (32)-(35) and (49), which cau be solved for the 
five unknowns: 7?i, <?i, K^, K i, and Mi. This is done numerically (see Appendix 2) and 
the results arc shown in the next section. 
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4. RESULTS 






4.1 Position of the Equilibrium Point 

As it can be seen from Figure 6 the equilibrium point is always ahead of the leading 
edge of the plate and it goes farther as the angle of attack increases reaching a maximum 
distance of approximately 53 percent of the chord from the leading edge (measured along 
the chord line) at 45 degrees angle of attack. 

At the same time (see Figure 7) after reaching a maximum distance above the chord 
line of approximately 16 percent of the chord (measured perpendicular to the chord line), 
it starts moving downwards crossing the chord line at a a 35° and getting farther under 
it at higher angles of attack. 


4.2 Vortex and Sink Strengths 

The bound vortex strength increases almost linearly with angle of attack while the 
leading edge vortex strength grows nonlincarly reaching a maximum at a = 45° and 
dropping to zero at a ~ 73° where the lower stagnation point reaches the trailing edge and 
moves olT the plate (Figure 8). Thus the total circulation is also nonlinear with angle of 
attack and has a maximum at a = 50° (Figure 10). The sink strength on the other hand, 
starting with very small values increases almost linearly with angle of attack (Figure 9). 

4.3 Pressure Distribution on the Plate 

The pressure distribution has been calculated at the plane (the most convenient 
one since the plate is horizontal and lies along the £'-axis). 

From the sequence of transformations shown in Figure 2 we have 


= + 


= ze ,a -f 


.-•a 


= re t0 e tn + — t 

r 


\o 


or 

a 2 

= r [cos(0 + a) -f- i sin(0 + a)] -|- — [cos(0 + a) — t sin(0 + a)]. 
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and since 


$' = £'+*V 

we liave 

cos(0 -f a) 
and 

rf' = - y) sin(^ + or) 

on the surface of the plate r = a so 

Cp = 2acos(0 + ar) 

«!» = ® 


(50) 

(51) 

(52) 

(53) 

(54) 


To find the velocity field around the plate, a similar procedure as in Section 2.3 yields: 

dw _ dw dz 
d? ~ !z'd? 

dw z* 

dz z 2 e iu — a 2 e~ ia 


or 


rV w 


K - *“».) = • («r. - <<*»,) • 

and on the surface (r = u) we get finally 


from which 


1 

COS0 

(55) 

Uf ‘» “ ~ 2 

sin(0 -1- a) U °* a 

1 

sin0 

(56) 

= ~2 

’ sin(0 -f at) U ° M 

w ap = ^Utf, e csc(0-f«) 

(57) 


Now from equation (5) and using equations (27)— (31) we find that 


= -t+ t -K+ aft * y ;;‘) + y?r 11 

Voo 1 + fif - 2i?i cos(0 - ) 


(58) 
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if m i c 1 )»«■**.> * £*■** x . 


J*_ 


;s. .. Jt ' i T 


Substituting equation (58) into (57) as well as the corresponding solution sets for .each 
angle of attack we can find the pressure distribution from 


^p.r — 1 - 


va 


(59) 


C v , p has been plotted for four angles of attack in Figure 11 where the horizontal axis 
measures the non-dimensional parameter 

£' 

n = cos(0 + a) = 

2 a 

which varies from —1 at the leading edge to +1 at the trailing edge. 

It can be seen from the plottings that both stagnation points move towards the 
trailing edge as the angle of attack increases but the lower one moves faster, being at 
about 70 percent of the chord from the leading edge at a — 60° while the upper one is 
only at the mid-chord point. 


4.4 Lift, Drag, and Pitching Moment 


At this point a comparison will be made between the results of: 

(i) the classical solution of totally attached (low over the plate with the Kutta condition 
satisfied at the trailing edge only and a singularity (infinite suction) at the leading 
edge (Figure 12a), 

(ii) this model with a detached vortex and a sink which make possible the satisfaction 
of the Kutta condition at both edges removing thus the leading edge singularity 
(Figure 12b), 

(iii) the Helmholtz solution of totally separated flow over the plate (Figure 12c), 

The lift, drag and moment coefficients for the three models arc given below: 


C Ll = 2 jt sin a 


Cl 2 = irK 


Clz 


it sin 2 a 
4 + 7T sin a 


Cdi — 0 Cmi = 7 sin 2a 

4 

C D 2 — -it Mi Cm 2 — t sin 2a 

4 

* 2n sin 2 a 
Cdz = j— — : — 

4 4- it sin a 


(00) 

( 01 ) 

( 02 ) 
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From Figure 13 it can be seen that our model with partially separated flow (as the 
leading edge vortex can be thought of), gives the highest lift coefficient at least up to 
of -60°. 

At higher angles it looks like the first model gives slightly higher lift coefficients, but 
this is misleading since at high angles of attack the flow separates at some point on the 
upper surface, resulting in significantly lower lift coefficients. Thus the first model breaks 
down at high angles of attack. 

From Figure 14 it can be seen that our model gives much higher drag coefficients 
even than the third model in which the flow is totally separated on the upper surface of the 
plate. This should not be surprising however, because the drag in our model is associated 
with a large momentum loss of all the fluid that dissappears into the sink (see Figure 5). 

As a confirmation to the above comes Figure IS which shows the same lift to drag 
ratio for the second and third models (since L/D = tan a in both cases). In other words 
since the second model yields a much higher lift than the third one, it also gives a much 
higher drag which can be thought of as induced drag (lift related drag). 

' Finally, it is interesting to note that our model has the same moment coefficient and 
the same position of aerodynamic center (at the quarter chord point) as the first model. 
The proof is as follows: Since there is a uniform wind, the velocity at a great distance from 
the plate must tend simply to the wind velocity, and therefore if \z\ is sufficiently large we 


may write 

dw ..A B 


T*= v ~ + 7 + ? 

or 

B 

tu = VqoZ + A In z + — 


Now the force and the moment on the plate can be found from the theorem of Blasius 
and it turns out (after performing the integration) that the force depends only on A while 
the moment depends only on B . However, comparing the complex potential for the first 
two models wc can sec that B is the same regardless of the presence of the vortex/sink 
combination and therefore the moment ought to be the same for the two models. In 
equations (GO) and (61) Cm is taken about the mid-point of the plate and is therefore 
positive. The same relations but with a minus sign on the right hand side arc valid for 
Cm about the leading edge. Equations (GO) and (Gl) arc plotted in Figure 16. 


5. CONCLUSIONS 


(i) The present inviscid, incompressible, two-dimensional model of a flat plate with a 

i 

detached vortex close to its leading edge indicates that lift coefficients up to around 
6 are achievable. Higher values should also be possible if thickness and camber are 
added, considering an airfoil instead of a flat plate (reference 3). 

(ii) In order to satisfy the Kutta condition at both the trailing edge and the leading 
edge, the presence of the sink is necessary (see equations (32) and (33)). This is 
in agreement with physical observations of the leading edge vortex which forms 
over delta wings at high angles of attack where the spanwise pressure gradient due 
to sweep angle evacuates the vortex core. It is also in agreement with the model 
presented in reference 3. 

(iii) For a given lift (Cl) there are two possible solutions for the location of the equi- 
librium point and the corresponding strengths of the vortex and sink. That of the 
lower angle of attack gives a weaker sink and therefore less drag, while that of the 
higher angle of attack gives a stronger sink and the associated higher drag. The vor- 
tex strengths do not differ much for the two solutions since they are closely related 
to the lift (which is the same for the two solutions). This result is also in agreement 
with reference 3. 

(iv) The upper limit found in the Cl versus a curve suggests that if a stronger vor- 
tex would exist at the equilibrium point, unrealistic supercirculation would occur 
resulting in the streamlines going entirely around the system. 

(v) A limitation of the present model appears at a esc 73°, above which the lower stag- 
nation point moves off the plate making thus imposc'blc the (low pattern depicted 
in Figure 5, on which this model is based. 

(vi) A comparison between Figures 8 and Ic shows very good agreement (at least qualita- 
tively) between the results derived here and those found experimentally in reference 
1 . 

(vii) Finally, the presence of the vortex and the sink docs not alTcct the position of the 
aerodynamic center which remains at the quarter chord point. 
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Figure 4 • Force diagram on the flat plate. 
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Figure 15. Lift to drag ratio versus angle of attack. 
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Figure 16 . Pitching moment characteristic. 




APPENDIX 1 


Al.l 

Consider a vortex of strength ki at the point zi outside of a cylinder \z\ — a. Then, 
the complex potential is . 


—/(*)+/ (7) 

= iki ln(z — zi) — iki In — zi^ 

according to the circle theorem (reference 2) w can also be written as 

w = iki hi z -f iki bi(* — zi) — iki In ^z — -f- const. 

which shows three vortices: one at the point Zi\ one at the point £ which is the inverse 
square point of Zi with respect to the cylinder; and one at the origin (center of the cylinder). 

A1.2 

Consider a sink of strength mj at the point Zi outside of a cylinder \z\ — a. Then, 
the complex potential is 


or 


w 


w 


w 


-/(•)+/( 7) 

(a? \ 

—■ mi ln(z — zi) + mi In l — — zij 
= —mi In z + mi ln(z — zi) -f mi In ^z — + const. 


equivalent to: a sink at the point zi, a sink at the point and a source at the origin. 


APPENDIX 2 


A four-dimensional Newton-Raphson algorithm has been used to solve the system 
of equations (32), (33), (34), and (35) for /Zi, Ky, Ki after elimination of Mi from 
equation (49). 


non 


ORIGINAL S3 
OF POOR QUAWO! 


PROGRAM NR4D 

REAL K0,K1,K,J2,J3,J4,J5,J6,J7,I1,I2,I3,I4,I5,I6,I7,I8,I9 


DATA 


ALFA=? 

ALF=ALFA* 3 . 14159/180 .0 
30 WRITE (5,110) 

110 FORMAT ( ' GIVE UNDERRELAXATION FACTOR') 

READ (5,120) C 
120 FORMAT (F) 

C T=THETA, R=R1 

R=? 

T=? 

K0=? 

Kl=? 

10 WRITE (5,100) R,T,K0,K1 

100 FORMAT (4E15.6) 

c 

T1=6*T+2*ALF 

T2=2*T-2*ALF 

T3=T+ALF 

K=K0+K1 

C— — 

A1=2*R**4* (R*R-1) * (9*R*R-5) *COS (T) *COS (Tl) 

A2=-4*R*R* (R*R-1) * (7*R*R-3) *COS (T) *COS (4*T) 

A3=2* (R‘R-1) * (5*R*R-1) *COS (T) *COS (T2) 

A4=-2*R**4* (9*R**4-5) *SIN(T) *SIN(T1) 

A5-4*R*R* (7*R* *4-3) *SIN(T) *SIN(4*T) -2* (5*R**4-1) *SIN(T) *SIN(T2) 
A6=-4*K0*R* *5* (4*R*R-3) *SIN(T1) +8*K0*R**3* (3*R*R-2) *SIN(4*T) 
A7=-4*K0*R* (2*R*R-1) *SIN(T2) +12*K1*R**5*SIN(T1) 

A8=2*K1*R* *3* (21*R*R-22) *SIN (4*T) +14*K1*R* *3* (3*R*R-2) *COS(4*T) 
A9=8*K1*R* (3*R*R-1) *SIN (T2) +12*K1*R* (2*R*R-1) *COS (T2) 

B1=6*R* (K1+K*TAN (ALF) ) * (COS (4*ALF) -SIN(4*ALF) ) 

B2=-12*R* *5*K*TAN (ALF) *COS (Tl) 

B3=3*K*R**3* (14*R*R-4) *TAN (ALF) *COS (4*T) 

B4=14*K*R* * 3* (3*R*R-2) *TAN (ALF) *SIN(T) 

B5=K*R*TAN (ALF) * (8* (3*R*R-2) *COS(T2) +12* (2*R*R-1) *SIN(T2) ) 

c 

All = A1+A2+A3+A4+A5+A6+A7+A8+A9+B1+B2+B3+B4+B5 

Q. 

B6=-2*R**5* (R*R-1) **2* (SIN (T) *COS (Tl) +6*COS (T) *SIN(T1) ) 

B7=4*R* *3* (R*R-1) * *2* (SIN(T') *COS (4*T) +4*C0S (T) *SIN(4*T)) 
B8=-2*R* (R*R-1) **2* (SIN(T) *COS (T2) +2*COS (T) *SIN(T2) ) 

B9=-2*R* *5* (R**4-l) * (COS (T) *SIN(T1) +6*SIN(T) *COS (Tl) ) 

Cl— 4*R**3* (R**4-l) * (COS (T) *SIN(4*T) +4*SIN(T) *C0S(4*T)) 

C2=-2*R* (R**4-l) * (COS (T) *SIN(T2) +2*SIN(T) *COS (T2) ) 
C3=-12*K0*R**6* (R*R-1) *COS (Tl) +16*K0*R**4* (R*R-1) *C0S(4*T) 
C4=-4*K0*R*R* (R*R-1) *COS (T2) +12*K1*R**6*C0S (Tl) 

C5=4*K1*R**4* (7*R*R-11) *COS (4*T) -28*K1*R**4* (R*R-1) *SIN(4*T) 
C6=4*K1*R*R* (3*R*R-2) *COS (T2) -12*K1*R*R* (R*R-1) *SIN(T2) 
C7=4*K*R* *4*TAN (ALF) * (3*R*R*SIN (T2) - (7*R*R-3) *SIN (4*T) ) 
C8=28*K*R**4* (R**2-l) *TAN(ALF) *C0S(4*T) 

C9=-4*K*R*R*TAN (ALF) * ( (3*R*R-4) *SIN(T2)-3* (R*R-1) *COS(T2) ) 

C 

A12 = B6+B7+B8+B9+C1+C2+C3+C4+C5+C6+C7+C8+C9 

C 

D1=-2*R**6* (R*R-1) *SIN^T1) +4*R**4* (R*R-1) *SIN(4*T) 

D2=-2*R*R* (R*R-1) *SIN(7’) 


G 

G 


G 

G 


C 

c- 


c 

c 


c 

c- 


D3=3* (R*R-1) *TAN (ALF) * (COS (4*ALF) -SIN(4*ALF) ) 
D4=-2*R**6*TAN(ALF) *COS (Tl)+R**4* (7*R*R-3) * TAN (ALF) *C0S(4*T) 
D5=7*R**4* (R*R-1) *TAN (ALF) *SIN(T) 

D6=2*R*R* (3*R*R-4) * TAN (ALF) *COS(T2) 

D7=6*R*R* (R*R-1) *TAN (ALF) *SIN(T2) 


A13=D1+D2+D3+D4+D5+D6+D7 


D8=2*R**6*SIN(T1)+R**4* (7*R*R-11) *SIN(4*T) 

D9=7*R**4* (R*R-1) *COS(4*T) +2*R*R* (3*R*R-2) *SIN(T2) 

E1=6*R*R* (R*R-1) *COS (T2) 

E2=3* (1+TAN (ALF) ) * (R*R-1) * (COS (4*ALF) -SIN (4*ALF) ) 
E3=-2*R**6*TAN (ALF) *COS (Tl) +R**4* (7*R*R-3) *TAN (ALF) *C0S(4*T) 
E4=7*R**4* (R*R-1) *TAN (ALF) *SIN(T) 

■E5=2*R*R*TAN(ALF) * ( (3*R*R-4) *COS (T2) +3* (R*R-1) *SIN(T2) ) 


A14=D8+D9+E1+E2+E3+E4+E5 


E6=2*R**4* (R*R-1) * (9*R*R-5) *COS(T) *SIN(T1) 

E7=-4*R*R* (R*R-1) * (7*R*R-3) *COS (T) *SIN(4*T) 

E8=2* (R*R-1) * (5*R*R-1) *COS(T) *SIN(T2) 

E9=2*R**4* (9*R**4-5) *SIN(T) *COS (Tl) 

G1=-4*R*R* (7*R**4-3) *SIN(T) *COS (4*T) +2* (5*R**4-1) *SIN(T) *COS (T2) 
G2=4*K0*R**5* (4*R*R-3) *COS (Tl) -8*K0*R**3* (3*R*R-2) *C0S(4*T) 
G3=4*K0*R* (2*R*R-1) *COS (T2) -12*K1*R* *5*COS (Tl) 

G4=-2*K1*R**3* ( (21*R*R-22) *COS (4*T) +7* (3*R*R-2) *SIN(4*T) ) 
G5=-4*K1*R* (2* (3*R*R-1) *COS (T2) +3* (2*R*R-1) *SIN(T2) ) 

G6=-6* (K1-K*TAN (ALF) ) *R* (COS (4*ALF) -SIN (4*ALF) ) 

G7=-12*K*R* *5*TAN (ALF) *SIN(Tl) 

G8=K*R * * 3 * TAN (ALF ) * (3* (14*R*R-4) *SIN(4*T) +14* (3*R*R-2) *COS (4*T) ) 
G9=4*K*R*TAN (ALF) * (2* (3*R*R-2) *SIN (T2) +3* (2*R*R-1) *COS (T2) ) 


A21 = E6+E7+E8+E9+G1+G2+G3+G4+G5+G6+G7+G8+G9 


I1=2*R**5* (R*R-1) **2* (-SIN (T) *SIN(T1) +6*C0S (T) *COS(Tl) ) 
I2=-4*R**3* (R*R-1) **2* (-SIN(T) *SIN(4*T) +4*C0S(T) *COS (4*T) ) 
I3=2*R* (R*R-1) * *2* (-SIN (T) *SIN (T2) +2*COS (T) *COS(T2)) 

14=2 *R* *5* (R* *4-1) * (COS (T) *COS (Tl) -6*SIN (T) *SIN (Tl) ) 

I5=-4*R* * 3* (R* *4-1) * (COS (T) *COS (4*T) -4*SIN (T) *SIN (4*T) ) 

16=2 *R* (R* *4-1) * (COS (T) *COS (T2) -2*SIN (T) *SIN (T2) ) 
I7=-4*K0*R*R* (R*R-1) * (3*R**4*SIN(T1) +SIN(T2) -4*R*R*SIN(4*T) ) 
I8=4*K1*R**4* (3*R*R*SIN (Tl) + (7*R*R-11) *SIN(4*T) ) 
I9=-4*K1*R*R* (7*R*R* (R*R-1) *COS (4*T) - (3*R*R-2) *SIN (T2) ) 
P1=-12*K1*R*R* (R*R-1) *COS (T2) 

P2=-4*K*R**4*TAN (ALF) * (3*R*R*COS (Tl) - (7*R*R-3) *COS(4*T)) 
P3=-28*K*R* *4* (R**2-l) *TAN (ALF) *SIN(4*T) 

P4=4*K*R*R*TAN (ALF) * ( (3*R*R-4) *COS (T2) -3* (R*R-1) *SIN (T2) ) 


A22 = I1+I2+I3+I4+I5+I6+I7+I8+I9+P1+P2+P3+P4 


P5=2 *R*R* (R*R-1) * (R**4*C0S (Tl) -2*R*R*COS (4*T) +COS (T2) ) 
P6=3* (R*R-1) *TAN (ALF) * (COS (4*ALF) -SIN (4*ALF) ) 
P7=-2*R**6*TAN (ALF) *SIN(T1) 

P8=R* *4*TAN (ALF) * ( (7*R*R-3) *SIN(4*T) +7* (R*R-1) *COS (4*T) ) 
P9=2*R*R*TAN (ALF) * ( (3*R*R-4) *SIN(T2) +3* (R*R-1) *COS (T2) ) 


A23=P5+P6+P7+P8+P9 


Ql=-R**4* (2*R*R*COS (Tl) + (7*R*R-11) *COS (4*T) +7* (R*R-1) *SIN (4*T) ) 
Q2=-2*R*R* ( (3*R*R-2) *COS (T2) +3* (R*R-1) *SIN(T2) ) 


c 

c 


c 


c 


c 

c 


c 

c 


c 


Q3=-3* (1-TAN (ALF) ) * (R*R-1) * (COS (4*ALF)-SIN(4*ALF) ) 
Q4=-2*R**6*TAN(ALF)*SIN(T1) 

Q5=R* *4*TAN (ALF) * ( (7*R*R-3) *SIN(4*T) +7* (R*R-1) *COS (4*T) ) 
Q6=2*R*R*TAN (ALF) * ( (3*R*R-4) *SIN(T2) +3* (R*R-1) *COS (T2) ) 


A24=Q1+Q2 +Q3+Q4+Q5+Q6 


A31=2* (2*SIN (ALF) -K) * (R-COS (T3) ) +2*K1*R+2*K*TAN(ALF) *SIN(T3) 
A32=2*R* (2*SIN (ALF) -K) *SIN(T3) +2*K*R*TAN(ALF) *COS(T3) 
A33=-(l+R*R-2*R*COS(T3) ) +2*R*TAN (ALF) *SIN(T3) 
A34=-(l+R*R-2*R*COS(T3))+R*R-l+2*R*TAN(ALF)*SIN(T3) 


A41=2* (2*SIN (ALF) +K) * (R+COS (T3) ) +2*K*TAN(ALF) *SIN(T3) -2*K1*R 
A42=2*K*R*TAN (ALF) *COS (T3)-2*R* (2*SIN(ALF) +K) *SIN(T3) 
A43=l+R**2+2*R*COS (T3) +2*R*TAN (ALF) *SIN(T3) 
A44=l+R**2+2*R*COS (T3) +2*R*TAN(ALF) *SIN(T3)-(R*R-1) 


FA=2*R* (R*R-1) **2*COS (T) * (R**4*C0S (Tl)-2*R*R*COS (4*T) +COS (T2)) 
FB=-2*R* (R**4-l) *SIN(T) * (R**4*SIN(T1) -2*R*R*SIN(4*T) +SIN(T2) ) 
FC=-2*K0*R*R* (R*R-1) * (R**4*SIN(T1) -2*R*R*SIN(4*T) +SIN(T2) ) 
FD=2*K1*R**6*SIN(T1) 

FE-K1*R**4* ( (7*R*R-11) *SIN(4*T) +7* (R*R-1) *COS(4*T) ) 
FG=2*K1*R*R* ( (3*R*R-2) *SIN (T2) +3* (R*R-1) *COS (T2) ) 

FH=3* (K1+K*TAN(ALF) ) * (R*R-1) * (COS (4*ALF) -SIN (4*ALF) ) 
FI--2*K*R**6*TAN (ALF) *COS (Tl) 

FJ=K*R**4*TAN (ALF) * ( (7*R*R-3) *COS (4*T) + 7 * (R*R-1) *SIN(4*T) ) 
FK-2*K*R*R*TAN(ALF) * ( (3*R*R-4) *COS (T2) +3* (R*R-1) *SIN(T2) ) 


FI = FA+FB+FC+FD+FE+FG+FH+FI +FJ+FK 


FL=2*R* (R*R-1) **2*COS (T) * (R**4*SIN(T1)-2*R*R*SIN(4*T) +SIN(T2) ) 
FM-2*R* (R* *4-1) *SIN (T) * (R**4*COS (Tl) -2*R*R*COS (4*T) +COS (T2) ) 
FN=2*K0*R*R* (R*R-1) * (R**4*COS(Tl)-2*R*R*COS(4*T) +COS (T2) ) 
FO=-2*Kl*R**6*COS (Tl) 

FP=-K1*R**4* ( (7*R*R-11) *COS (4*T) +7* (R*R-1) *SIN(4*T) ) 
FQ=-2*K1*R*R* ( (3*R*R-2) *COS (T2) +3* (R*R-1) *SIN(T2) ) 

FR=-3* (K1-K*TAN (ALF) ) * (R*R-1) * (COS (4*ALF) -SIN(4*ALF) ) 
FS=-2*K*R* *6*TAN (ALF) *SIN (Tl) 

FT=K*R**4*TAN (ALF) * ( (7*R*R-3) *SIN (4*T) +7* (R*R-1) *COS (4*T) ) 
FU=2*K*R*R*TAN (ALF) * ( (3*R*R-4) *SIN (T2) +3* (R*R-1) *COS (T2) ) 


F2 = FL+FM+FN+FO+FP+FQ+FR+FS+FT+FU 


FV=(2*SIN (ALF) -K) * (l+R*R-2*R*COS (T3) ) 
FW=K1* (R**2-l) +2*K*R*TAN(ALF) *SIN(T3) 
F3=FV+FW 

FX= (2*SIN (ALF) +K) * (l+R*R+2*R*COS (T3) ) 
FY=2*K*R*TAN (ALF) *SIN(T3)-K1* (R*R-1) 
F4=FX+FY 


S1=A11*A22* (A3 3 * A44-A43 * A34) +A11*A23* (A32*A44-A42*A34) 
S2=A11*A24* (A32*A43-A42*A33)+A12*A21* (A33*A44-A43*A34) 
S3=A12*A23* (A31*A44-A34*A41) +A12*A24* (A31*A43-A41*A33) 
S4=A13*A21* (A32*A44-A34*A42) +A13*A22* (A31*A44-A34*A41) 
S5=A13*A24* (A31*A42-A41*A32) +A14*A21* (A32*A43-A33*A42) 
S6=A14*A22' V (A3 1 * A43-A33 * A41 ) +A14*A23* (A31*A42-A41*A32) 
DET = S1+S2+S3+S4+S5+S6 


C 


J2=F1*A22* (A33 * A44-A43 * A34) +F1*A23* (A32*A44-A42*A34) 
J3=F1*A24* (A32 *A43-A42 *A33) +A12*F2* (A33*A44-A43*A34) 







J4=A12*A23* (F3*A44-A34*F4) +A12*A24* (F3*A43-F4*A33) 
J5=A13*F2* (A32*A44-A34*A42)+A13*A22* (F3*A44-A34*F4) 
J6=A13*A24* (F3*A42-F4*A32)+A14*F2* (A32*A43-A33*A42) 
J7=A14*A22* (F3*A43-A33*F4) +A14*A23* (F3*A42-F4*A32) 
HI = (J2+J3+J4+J5+J6+J7) / DET 


U1=A11*F2* (A3 3 * A44-A4 3 * A34) +A11*A23* (F3*A44-F4*A34) 
U2=A11*A24* (F3*A43-F4*A33) +F1*A21* (A33*A44-A43*A34) 
U3=F1*A23* (A31*A44-A34*A41) +F1*A24* (A31*A43-A41*A33) 
U4=A13*A21* (F3*A44-A34*F4) +A13*F2* (A31*A44-A34*A41) 
U5=A13*A24* (F4*A31-F3*A41) +A14*A21* (F3*A43-F4*A42) 
U6=A14*F2* (A31 *A43-A33*A41) +A14*A23* (F4*A31-F3*A41) 
H2 = (U1+U2+U3+U4+U5+U6) / DET 


V1=A11*A22* (F3*A44-F4*A34) +A11*F2* (A32*A44-A42*A34) 
V2=A11*A24* (A32*F4-A42*F3) +A12*A21* (F3*A44-F4*A34) 
V3=A12*F2* (A31*A44-A34*A41)+A12*A24* (A31*F4-A41*F3) 
V4=F1*A21* (A32*A44-A34*A42) +F1*A22* (A31*A44-A34*A41) 
V5=F1*A24* (A31*A42-A41*A32) +A14*A21* (A32*F4-F3*A42) 
V6=A14*A22* (A31*F4-F3*A41) +A14*F2* (A31*A42-A41*A32) 


H3= (V1+V2+V3+V4+V5+V6) / DET 

C | 

W1=A11*A22* (A33*F4-A43*F3) +A11*A23* (A32*F4-A42*F3) | 

W2=A1*F2* (A32*A43-A42*A33) +A12*A21* (A33*F4-A43*F3) 1 


W3=A12*A23* (A31*F4-A41*F3) +A12*F2* (A31*A43-A41*A33) 
W4=A13*A21* (A32*F4-A42*F3) +A13*A22* (A31*F4-A41*F3) 
W5=Al3*F2* (A31 * A42-A41 *A32) +F1*A21* (A32*A43-A33*A42) 
W6=F1*A22* (A31*A43-A33*A41) +F1*A23* (A31*A42-A41*A32) 
H4 = (W1+W2+W3+W4+W5+W6) / DET 


R =R-C*H1 
T =T-C*H2 
K0=K0-C*H3 
K1=K1-C*H4 

IF (ABS (HI) .LT. 1 . E-3.AND.ABS (H2) .LT. 1 . E-3) GO TO 20 
IF (ABS (H3) . LT . 1 . E-3 . AND . ABS (H4) . LT . 1 . E-3) GOTO 20 
GO TO 10 

20 WRITE (5,100) R,T,K0,K1 
WRITE (5,130) 

130 FORMAT ( * IF YOU WANT ANOTHER RELAXATION FACTOR, TYPE 1*) 
READ (5,140) IFF 
140 FORMAT (I) 

IF (IFF.EQ.l) GO TO 30 

STOP 

END 



